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Abstract 

An analytical approach to heat transfer for crystal growth 
in a Bridgman-Stockbarger configuration has been developed. All 
first order effects on the axial temperature distribution in a 
solidifying charge are analyzed on the basis of a one-dimensional 
model whose validity could be verified through comparison with 
published finite difference analyses of two-dimensional models. 

The present model includes an insulated region between axially 
aligned heat pipes and considers the effects of charge diameter, 
charge motion, thickness and thermal conductivity of a confining 
crucible, thermal conductivity change at the crystal-melt inter- 
face, generation of latent heat at the interface and non-infinite 
charge length. Results are primarily given in analytical, form and 
can be used without recourse to computer work for both improved 
furnace design and optimization of growth conditions in a given 
thermal configuration. 
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1, INTRODUCTION 


Conventional technique's for crystal growth from the melt 
prove increasingly inadequate in meeting property requirements 
for electronic materials as dictated by advanced device technology. 
In this context it is of interest that renewed attention is 
recently being paid to seeded vertical Bridgman growth where the 
driving force for interfering free melt convection is reduced 
and the critical axial and radial thermal gradients are at least 
in principle readily controllable. Motivation for a 
thorough re-examination of this crystal growth technique is 
provided by ample evidence in the open literature that basic heat 
transfer considerations have been largely ignored in its applica- 
tion. These considerations are now recognized to strongly affect 
both crystalline and chemical perfection of the resulting solid. 

In recent years several thermal analyses of Bridgman-type 

crystal growth systems have been reported in the literature. 

Noteworthy among these is a series of outstanding publications by 

Wilcox et al. ' ' ' in which the effects of several dimensionless 

parameters on both axial and radial temperature distribution were 

analyzed and significant conclusions were drawn. The modeled 

systems, however, were simplified in several respects and the two- 

dimensional formulation (except for part of Chang and Wilcox 

precluded simple analytical results. Of interest is also the work 
. (4) (5 6) 

of Davis ' and dyne ' ' which demonstrates that one-dimensional 
models can accurately predict experimental axial temperature vari- 
ations. Their models were not nondimensional and, consequently, 
significant thermal and geometric parameters which control the 


system behavior were not identified, and their conclusions cannot 
readily bo extended to systems of different parameter values. 

The present work is part of a comprehensive integrated theo- 
retical and experimental approach to vortical Bridgman growth in 
which one- and two-dimensional heat transfer models are used as a 
basis for optimized system design. The system ‘(presently in con- 
struction) makes use of high temperature heat pipes which provide 
for axial and axi-symmetric temperature uniformity and thus for 
thermal boundary conditions which permit a meaningful theoretical 
analysis. The configuration (Fig. 1) includes, between aligned 
heat pipes, a gradient control region in which the crystal-melt 

interface is to be located. The purpose of this region, first 

( 1 ) 

suggested by Chang and Wilcox, is to provide control of heat 
•flow near the interface, required for the establishment of desired 
radial thermal gradients. The design and operation parameters are 
expected to provide for the system under construction a wide range 
of critical axial gradients and radial gradients leading to growth 
interface morphologies ranging from convex to concave. Primary 
attention is focused, hov/ever, on the establishment of growth 
conditions which preclude nucleation at the confining boundary 
with the simultaneous minimization of radial segregation effects. 
The use of a magnetic field for the reduction of free melt convec- 
tion is provided for. 

This study analyzes the effects of both system and operation 
parameters on the axial temperature distribution in the charge on 
the basis of a one-dimensional model.* This particular approach 

*Radial temperature variations in the charge are the subject of a 
publication (Part II) presently in preparation. 
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was selected since oversimplified models are not required in order 

to provide useful and easily applied analytical results. The heat 

transfer analysis focuses on the following effects; geometry and 

material of the charge-confining crucible, diameter effect of the 

charge, length of the gradient control region, charge motion, 

thermal coupling to the furnace, charge 'length, generation of 

* 

latent heat at the crystal-melt interface, and thermal conductivity 
change at the interface. Special attention is given to the axial 
temperature gradient at the interface and the axial position of the 
interface within the gradient control region. 

The purpose of this analysis, which considers all first order 
effects on axial thermal gradients, is to provide a basis for the 
design of Bridgman-type growth systems and to give guidelines for 
optimized axial gradient control during the execution of growth 
experiments. Its primary asset is the timesaving element, sines 
meaningful information can be obtained from the analytical results 
without resort to computer work. 

2. DEVELOPMENT OF HEAT TRANSFER MODEL 

The factors of concern for the development of a one-dimensional 
heat transfer model of the Bridgman growth system depicted in 
Fig. 1 are shown in Fig. 2. Hot and cold heat pipes comprise the 
hot and cold zones; the region between them is called the gradient 
zone. The length of the charge is broken down into L^j, L^, and L^ 
within the hot, gradient and cold zones respectively. The charge 
is lowered through the furnace with a velocity V, has liquid and 
solid portions with different thermal conductivity, and has a 
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crystal-molt interface which generates latent heat. A crucible 
provides containment for the charge, 

'fhe model makes the following assumptions: 

1. Hot and cold zone furnace temperatures are 
uniform, reflecting the heat pipe action, 
and extend to infinity in cither direction. 

^ V * 

2. The system is at all times in a quasi-steady 
state; i.e. transients are neglected. 

3. No heat transfer by convection in the melt, 

4 . Heat exchange between the furnace and the charge 
is described by a heat transfer coefficient, h, 
that is constant within each zone; for simpli- 
fied models in which the crucible is not taken 
into consideration, h is calculated between the 
furnace and the surface of the charge; when a 
crucible is included in the modeling, h is calcu- 
lated between the furnace and the outer crucible 
surface. The gradient zone is assumed adiabatic, 
i.e. hg = 0. 

5. The one-dimensional model considers only axial 
heat transfer within the charge; i.e. the temper- 
ature in the charge is not a function of radius. 

With these assumptions, the charge is analagous to the moving 

( 7 ) 

thin rod treated by Carslow and Jaeger and the equation, in non 
dimensional form, describing the axial temperature distribution is 

2 

^ - Pe - 4Bi0 = -4Bi0^ tl] 


where : 

^ - z/D (non-dimensional axial coordinate measured 
from the center of the gradient zone) 
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0 B (non"<aimensional temperature) 

Of - non-dimensional furnace temperature 

Pe « Peclet number = VD/a 
Bi = Biot number = hD/k 

I 

a = thermal diffusivity 
h “ heat transfer coefficient 
k = thermal conductivity 

V = charge displacement rate (lowering rate) 

D = charge diameter 

Equation [1] is applied to each region of uniform properties and 
Bi. In the present model there are four such regions when placing 
the growth interface in the gradient zone: the hot and cold zones 

and the liquid and solid parts of the gradient zone. Boundary con- 
ditions are equality of temperature and continuity of flux between 
adjacent regions. At the crystal-melt interface, the flux condition 
gives: 


where 


^ liquid solid H dz 




^liquid^^solid 

, solid^^ref 


AH 


SL 


latent heat of solidification 


"p, solid 


= specific heat of the solid 


[ 2 ] 


(The Peclet number in eq. [2] is correctly based on the actual 
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growUh ratoond it is assumed here that the growth rate is equal to 
the lowering rate V.) 

For each region, eq. [1] yields two cKponential terms for the 
homogeneous solution and a particular solution that depends on 0^;. 
The constants of integration on the exponential terms are found 
using the boundary conditions. This has^ been done analytically 
for systems of infinite charge length; expressions for the axial 
gradient in the liquid and for the temperature at the crystal-melt 
interface are given in the Appendix. (If the complete axial temper- 
ature distribution is required, it is more convenient to determine 
the constants of integration by computer.) 

In consideration of the complex nature of the thermal model, 
a system with idealized parameter values, presently called a 
"symmetric" system, is used as a reference against which the effects 
of individual parameters will be assessed. A symmetric system is 
defined to have the following parameter values: 

1. Equal Biot number in hot and cold zone; i.e., 

B r H “ ^ c " ® ^ * 

2. Pe = 0. 

3. Equal solid and liquid thermal conductivities 
(Rl< = 1) • 

4. No generation of latent heat (R^j = 0). 

5. Equal charge lengths in the hot and cold zones 

(Ljj ^ L^) . 

Under such conditions, a symmetric system will have the 
axial temperature distribution, 


0(“r;) t Q(C) 


1 


[3] 
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3.1 Crucible Effect 


In Bridgman growth, all charges are confined in a crucible 
which, depending on conditions, varies in dimension and composition, 
containment of the charge tends to decrease axial gradients (Sen 
and Wilcox^^^). A crucible of low thermal conductivity lowers the 
gradient by adding thermal resistance between the charge and fur- 
nace, effectively decreasing Bi; one of high thermal conductivity 
lowers the gradient by "short-circuiting" heat flow within the 
crucible itself. 

The domain of eg. [1] is here considered to be the charge only 
and the crucible is not explicitly included. However, a simple 
model of the heat transfer within the crucible allows for its 
consideration in eg, [1] through modified Biot and Peciet numbers. 

The temperature distribution in the crucible obeys the heat 
conduction eguation in cylindrical coordinates. Neglecting Pe 
effects, we obtain; 


1 

"p 


.L 

3p 


(P 


30 


cr 


9p 


) + 


9 0 


cr 


= 0 


95' 


[4] 


where ; 

©cr = non-dimensional crucible temperature 

p = non-dimensional radial coordinate, r/D 
r = radial coordinate 

The first term of eg. [4] accounts for the radial thermal resis- 
tance of the crucible, while the second accounts for the "short- 


circuit". If, at each axial location of the crucible, 3^0 /9C^ 

cr ^ 

is considered to bo independent of p and equal to d^O/d?^ of the 
charge, eg. [4] can be integrated to yield a radial distribution 
of 0 . The boundary conditions are; 


do 

“‘•Of - Ocr> “ >^or -df 

♦•(at the outer cru- 
cible surface) 

(5a] 

« 0 (charge) 

(at the inner cru- 
cible surface 

[5b] 


where ; 

her “ heat transfer coefficient betV'/een furnace and 
outer curcible surface 

k « thermal conductivity of the crucible 
cr 

(Equation [5b] assumes that the crucible and charge are in intimate 
contact.) The heat flow from the crucible to the charge can thus 
be obtained in terms of d0 /dC and used to reformulate eg. [1] 
through "effective" Biot and Peclet numbers: 

Bi 

Bi = T n 1 ^ 5 

1+Bi -~in(S] +K(6 -1) tj~Bi Anfi 


Pe 


eff 


Pe[l +||Bi^^£n6] 

1 + Bi^j. [j{ 6^-1) - |in6] + K (6^-1) + 1 1 


[6b] 


®^eff 

^®eff 


« Biot to be used in eg. [1] 

=- Peclet to be used in eg. [1] 

= outside diameter of the crucible 


where : 


k * conductivity of the charge 
or liquid) 

6 - 


K - 


(either solid 


The relationship between the effective Biot number and the 

conductivity ratio K is shown in Fig. 3 for 6 * 1.25 and various 

Bij^j.. It can be seen that Bi^^^ is significantly reduced by both 

low and high values for K. (The relationship between Bi and axial- 

gradients will be established in Section 3.3; it shows that the 

present findings are in basic agreement with those of Sen and 
( 2 ) 

Wilcox. ' ) The conductivity ratio K which maximizes Bi^^^ can be 
obtained from eq. (6] j 


K 


(max Bi) 


Bi'^In^ ■ 
- 1 ) 


[71 


Taking typical values of Bi^^j. and 6, the conductivity ratio pro- 
viding for maximum axial gradients is found to lie bewteen 0.1 and 
1.0. (Since all charges require confinement, the Biot and peclet 
numbers used in subsequent sections are to be considered as Bi^^^ 
and unless otherwise stated.) 


3.2 Peclet Number Effect 


Pe is the ratio of axial heat transfer in the charge due to 
motion and due to conduction. The effect of Pe on the axial temper- 
ature distribution in an otherwise symmetric system, obtained 
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through eq. [X], is shown in Pig. 4, Xargt (o.g» high Xowering 
rate or Xow thermaX conductivity of the. charge) increases the tem- 
perature of the charge everywhere, As a consequence, in systems 
with constant moXting point temperature , the cryataX-meXt inter- 
face wiXX move toward the coXd zone? aXtornativoXy , the furnace 
temperuturos nwst be Xowered if the interface position is to 
remain fixed. 

In conventionaX growth exporiments, the Po effect is smaXX 
since the Xowering rates are smaXX (o.g. O.X-XO ym/sec) . A test 
criterion for its roXativc magnitude is provided through the 
characteristic roots of eq. [X] : 


m. 


•[X * (X + 



It can be seen that Pe disappears from eq. [8] when 


[81 


Po 


>> X, 


[93 


Equation [91 agrees with Chang and WiXcox^^^ who reported that the 
Po effect was stronger for smaller Bi, If the inequality in eq. 

[93 holds, the Pe effect is smaXX enough to satisfy the PoeXet 
number cirterion for symmetric systems (i.e. Po « 0). It should 
be noted (eq. [21 ) that the generation of latent heat at the growth 
interface is also dependent on Pe. Thus, Pe has a smaXX effect on 
the axial temperature distribution only if eg. [91 holds and the 
latent heat effect is small (see Section 3,6). 
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3.3 Blot Number Effisct 

The Biot number I through the heat transfer coefficient h, is 
a direct measure of the thermal coupling between the charge and 
the furnace. Typical values for the effective Biot number (i.e. 
after accounting for the crucible) vary from 0.05 for high conduc- 
tivity materials such as Ge to 5,0 for low conductivity materials 
such as CdTe, Axial temperature orofiles for several Bi, calcu- 
lated from eq. [1] for symmetric systems, are shown in Pig. 5. in 

n \ 

agreement with Chang and Wilcox, ' ' it is found that the charge 
temperature follows raore closely the furnace temperature and, as 
a result, the axial temperature gradient in the gradient zone 
increases as Bi increases. 

The axial gradient behavior at the crystal-melt interface 
presented in the Appendix can be simplified for symmetric 
systems: 


dO „ -1 


[ 10 ] 


where : 

X « non-dimensional zone length, L/D 
Bi = Bi„ “ Bi^ for symmetric systems 

n 

Xh = X^^ * infinity 

The dependence of the gradient (dO/di^) on Bi, according to eq. [10], 
is plotted in Pig. 6 for various gradient zone lengths (Xq) . The 
graphs show that the dependence of the axial gradient bn the gradi- 
ent zone length becomes stronger with increasing Bi. It can also 


t 


be seen (cu^ve Xq w 0) that there exist* a minimum Bi for any 
desired rton-dimonsional axial gradient. An obvious alternative 
to increasing 'Bi for achieving larger axial gradients is to in- 
crease the temperature difference between the hot and cold zones 
since: 


I 


dT 


A rrt 

3 ; 


lU) 


This approach r however, is contingent on the absence of adverse 
side effects such ns the development of excessive vapor pressures 
associated with increases of the hot zone temperature; it is also 
no viable alternative if, as a consequence of furnace temperature 
changes , the non-dimensional interface temperature is altered and 
results in the interface being shifted into an undesirable region 
within the gradient zone. 

The Biot numbers of the hot and cold zones are generally not 
equal due, for example, to the temperature dependence of radiative 
heat transfer. Per such conditions the zone with the larger Bi 
will more strongly influence the overall temperature level (Pig. 
7). Compared to the symmetric case (Bij^ ® Bi^^) , the temperature 
increases everywhere within the charge when Bij^ > Bi^^ and 
decreases everywhere within the charge when Bi„ > Bi„, In Pig, 7, 

w ii 

the crystal-melt interface location for each curve is placed in 


the center of the gradient zone and the unequal Bi effect appears 
as a change in nondimensional interface temperature. To retain 
the interface in this location during a growth experiment necessi- 
tates a lowering of the hot and/or cold zone temperature as Bin 
decreases. 
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The expression for the axial temperature gradient in the gra- 
dient zone (analagous to eg. [10 J for Bi^ Bij^) isv 

a? 

where ; 

Xr = “ infinity 

System is symmetric except Bi^ 7 ^ Bi^ 

Interface is at the center of the gradient zone. 

Equation [12] indicates that if one of the Biot numbers is much 
smaller than the other, it will control the axial gradient. Effort 
to increase axial gradients by adjusting Bi should therefore first 
be directed at the zone with the lower Bi. 

3.4 Charge Diameter Effect 

The effect of the charge diameter D on the axial temperature 
gradient is assessed on the basis of eg. [10] . In practice, chang- 
ing D will affect other parameters as well; a corresponding change 
in the geometry of the furnace cavity, for example, will modify 
both conduction and radiation heat transfer coefficients; changes 
in and 6 will also alter the effective Biot number. In order 
to permit isolation of the diameter effect, all other parameters 
are assumed fixed in the present analysis. 

Since the axial coordinate in eg. [10], is non-dimensionali 
with respect to D, it is considered more informative to compare 
axial gradients based on the dimensional axial coordinate z. If 
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tho gradient zone length, remains unchanged (i.e., varies 
as D”^) , the limiting cases for large and small Bi are given by: 


«L rLm 

Bi ->• « 
Xq i- 0 


(d0/dz)2 

(do7dz)T 


« 1 


[13a] 


lim 
Bi ■> 0 


“(de/dz) 2 ~ 

s 

f°ll 

(dO/dz) 


[°2l 


1/2 


[13 b] 


where : 

subscript = 1 denotes value before a change in diameter 
subscript « 2 denotes value after a change in diameter 

If Xq remains unchanged L« varies as D) , then eq. [10] gives 

VUV/ U/i ) T 


lim 

Bi ’>■ “ 
Xg ^ 0 




15: 


[14a] 


lim 
Bi ->• 0 


'(dO/dz) 2 
"(W/dz ) , 


D. 


1/2 


[14b] 


The relationship between the axial gradient and changes in D as a 
function of the initial Biot number, Bi^^, (i.e., Bi based on the 
initial system with charge diameter D^) , is presented in Pig. 8. 
Accordingly, increasing the charge diameter (D) will decrease the 
dimensional axial gradient unless compensated for by other changes. 
Tho effect of decreasing the charge diameter, on the other hand, is 
a function of Bi. For large Bi a decrease in D will produce a 
largo increase in axial gradient if associated with a simultaneous 


decrease in 


G' 


If Bi is small, regardless of the magnitude of X 


any decrease in D will enhance the axial gradient*, however, to a 
lessor extent than with a large Di. 

3.5 Thermal Conductiv ity Ch ange at the Crystal-Holt Intor^acG 
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The effect of on the axial temperature distribution according 

« 

to eq. [1] is shown in Fig. 9 for systems which arc otherwise 
symmetric (i.e. , Bi^j « Bi^/R^^) . In all instances, the charge phase 
with higher thermal conductivity exhibits a lower axial gradient 
because of lower thermal resistance to heat transfer in the axial 
direction. For charges with Rj^ > 1, the axial gradient in the 
molt near the interface is therefore less than in the solid. The 
functional relationship between the axial gradient in the melt at 
the mo3 t~crystal interface and Rj^, analogous to eg. [10] , is; 

<If> “ 

in melt at Bi/^ (R,, t rJ/ ) + X„(R„ + 1) 
interface ^ 


where ; 


Bi 


H 




R 


K 


B r 

« « infinity 

5 >f 1; system is otherwise syiranetric 


Interface is at the center of the gradient zone 


Comparing cqs. [15] and [10], it is seen that, independent of Bi, 
any increase in Rj^ results in a decrease of the axial gradient in 
the melt. The axial gradient as a function of R^^ according to 
eq. [15] is plotted in Fig, 10 for several Bi^ values with Xq " 1.0. 
Applying eq. [IZl to germanium (R^^ « 2.5), for example, it is found 
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that tho conductivity effect reduces the intorfacial melt gradient 
by about 50%. 


3.6 Latent Heat Effect 


The generation of latent heat of solidification at the melt- 

V 

crystal interface enters the solution of eq, [1] only through tho 
boundary condition at this interface (eq, [2]). (Pe in this section 
is not PQgjfff hut should be interpreted as defined in eq. (1].) 

For the symmetric case (Rj^ = 1) , the axial temperature gradient 
in tho liquid at the crystal-melt interface thus becomes i 


d0 

d? 


sr? 


-2 


H- Pe^^Rjj(Bi“^/^ + Xq) 
2 (Bi“^^^ + Xp) 


[16] 


i.vhere ; 

Bi - Bi^ “ Bijj 
Xj. “ " infinity 

PegRjj 7 ^ 0; system is otherwise synunetric 
Interface is at the center of the gradient zone 


For R,, 7 ^ 1, the axial gradient assumes the form: 

i\ 



Bi 


_-2 + (61 -^/=’- t X^) 


[16a] 


B^H ” 


where : 
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From eq. tl6] is can bo soon bhat the effect of the latent heat on 
axial gradients (and also on the temperature profile) will be small 
if : 


% 


^(Pe 


R„) (Bi"^/^ + \ q ) 5 <|) << 1 


S^H 


[17] 


Equation [17] demonstrates that the effect of latent heat on the 
axial temperature behavior increases as Bi decreases. 

Axial temperature profiles calculated from eq. [1] for various 
values of ^ are plotted in Pig. 11 where, in order to isolate the , 
latent heat effect, Pe was chosen small enough so that eq. [9] is 
satisfied. It can be seen that the generation of latent heat 


both increasas tho charge t<;3mperatuirt^ and decireases the axial gra- 
dient at the interface in the molt. For small values of (j>, the 
latent heat effect disappears and the axial temperature profile 
approaches that of symmetric systems. 


3 . 7 *' Infinite" Charge Length 

The contribution from the posit.' '•e characteristic root, eq. [8], 
to the solution of eq. [1] is normally small and is zero for an 
infinite charge length. The charge thus appears infinite in length 
when the contribution from the negative root also becomes small, 
i.e., for largo 5 , The temperature change within the hot or cold 
zone reaches approximately 99% of its final value when exp[(m_)5] ~ 
0.01, Using this as a criterion for infinite length; 



^ m 


[18] 
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whore : 

* length of charge in hot or cold zone for charge 
to appear infinitely long 

If Pe is small so that eg. [9] is satisfied^ the characteristic 
roots become: 

« ±2/Bi ' ‘ [19] 

and by substitution wo find: 



which is a useful expression for determining 

For any charge length which appears infinite in both the hot 
ana cold zones, tno temperature rield of the charge will not change 
(interface position remains fixed) as it is lowered through the gra- 
dient zone. Figure 12 shows the progression of axial temperature 
profiles as the charge moves from the hot zone to the cold zone for 
charge lengths less than infinite. It can bo seen that charge 
temperatures are displaced toward the hot zone furnace temperature 
when most of the charge is in the hot zone and vice-versa. Accord- 
ingly, to achieve constant interface position for non-infinite 
charge lengths, the nondimensional melt tompex-aturo must be reduced 
as the experiment proceeds. At constant nondimensional melt tem- 
perature, the interface will move from the cold zone toward the hot 
zone and the interface growth velocity will bo greater than the 
lowering rate. This finding is in _ agreement with the experimental 
results of Clyne^"’'^^ which suggest that longer charges and higher 


Bi will bond to stabilize the interface position (growth rate is 
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the same as the displacement rate) . 

When the charge length is not "infinite", appropriate boundary 
conditions must be applied to the ends of the charge. For a solid 
pull rod, an approximate Bi can, for example, be obtained by 
treating it as a simple fin, exchanging heat with the environment. 
The curves of Fig. 12 were calculated uaing for the ends the same 
Bi as for the circumference. 


3 . 8 Effect of Radial Gradients on the Axial Temperature Profile 

In the development of eq. [1] , the radial temperature distri- 
bution of 0 at each axial location ^ was considered constant. 
Presently, radial gradients in the charge are considered insofar as 
they affect the preceding results for the axial temperature distri- 
bution . 

The heat transfer between the furnace and charge accounted for 
in the Bi term of eq. [1] depends on the temperature difference be- 
tween the furnace wall and the charge surface. This temperature 
difference can be accommodated in eq. [1] ) 




-4Bi(0^ - Gg) 


[ 21 ] 


where : 

0g = surface temperature of charge 

Integrating each term of eq. [21] with respect to the cross sectional 
area of the charge and assuming that the average of the derivatives 
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is approximately equal to the derivative of the averages, we find? 


^ - Pe if « -4Bi[(0« - 0) 

dr ^ 


(0g"O)] 


[ 22 ] 


where ; • *• 

t 

X ^ ^ ’ <i(area) 

'area 


[In eq. [22] the radial gradients are expressed by the term (0 -F),] 

s 

If the radial temperature distribution within the charge is approxi- 
mated as was done for the crucible in Section 3.1, we obtain: 


(0 


s 


- Q) 


32 


[23] 


Combining eq. [23] with eq. [22] leads to a ‘relationship for an 
effective Bi, r which accounts for radial temperature gradi 

ents within the charge : 


Bi 


* 

ef f 


Bi 


1 + 


Bi 


[24] 


where: 


Bi = Biot as defined in Section 3.1 [accounting 
for crucible effect, but not for radial 
gradient effect) 

Bi*^^ = Biot accounting for both crucible effect and 
radial temperature gradients in the charge 


Equation [24] can be used to approximate the effect of radial gradi- 
ents on the axial temperature distribution if Bi/8 is not small 
compared tu unity. If Bi/8 is small compared to unity, on the other 
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hand, eq. [1] reaj«onabXy predicts the axiaX temperature prefix© , 

The vaXidity of eq, [24] can be tested on the basis of data 
obtained through the two-dimensionaX finite difference modeX by 
Fu and WiXcox. fpiqure 3 of their paper presents the radiaX 
variation of the axiaX gradient at the interface of symmetric sys- 
tems for severaX Bi and X^, Point vaXues from these curves were 
numericaXXy integrated to obtain the average axiaX gradient over 
the cross section and were compared to the vaXues obtained using 
eqs. [XO] and [24] (see TabXe X). Considering the approximations 
required in the derivation of eg. [24] , the agreement of the data 
must be taken as excoXXent. The comparison demonstrates that the 
onc-diraensionaX modeXs presentXy used provide a moaningfuX repre- 
sentation of the parameters governing the axiaX temperature behavior 
of charges in a Bridgman-type growth configuration. 

4. DISCUSSION AND CONCLUSIONS 

The one-diraensionaX model developed in Section 2 has been shown 
to correlate weXX with a corresponding two-dimensionaX modeX of 
Sen and WiXcox concerning the axial temperature distribution in 
a solidifying charge. It is expected, therefore, that the results 
derived from this model accurately demonstrate the effects of 
furnace and material property parameters and can be applied to the 
optimization of both furnace design and execution of growth experi- 
ments. The model is easily applied to machine computation. Its 
primary asset, however, lies in the fact that it yields extremely 
useful and relatively simple analytical relationships. For example, 
criteria could thus be derived which define the conditions under v-/hich 
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the effects of charge motion, length of charge in the hot an(3 cold 
furnaces, and the generation of latent heat can be neglected; the 
effects of thickness and thermal conductivity of a confining cru- 
cible could be accounted for by a simple modification of the Bi 
parameter* . 

The model indicates that the axial .gradient in the liquid at 
the interface is adversely affected (decreases) by the following 
charge properties; 

•Large thermal conductivity 

•Large latent heat of solidification 

•Largo O'uquia/'^ solid’ 

Maximization of axial gradients in the charge can be approached 
through several operational and design options, each of which, how“ 
ever, has its limitations and potential drawbacks; 

Biot Number ; Axial thermal gradients can be increased by increasing 
the heat transfer coefficients (h) in both the hot and the cold 
zones (eq. [10]). Since in typical high temperature growth experi- 
ments heat transfer is largely controlled by radiation, it is imper- 
ative that furnace emissivities be kept high. (Within the gradient 
zone, however, undesirable radial heat transfer can be reduced 
through the installation of highly reflecting radiation shields.) 
Given the third power temperature dependence of the radiation heat 
transfer coefficient, an increase of the furnace temperature will 
also increase Bi. This approach, however, has its limitations in 
systems which develop high vapor pressures and thus require special 
crucible materials and construction, v;hich in turn may reduce Bi. 


24 


Heat conduction across the gap betv/oen the furnace and the cru- 
cible can significantly contribute to the overall heat transfer coef 
ficient if the gap width is small and the gas has a high thermal con 
ductivity. For example, He in a one millimeter gap will transfer 
about the same quantity of heat by conduction as is transferred by 
radiation at 900 ®C. Small gap widths, however,, accentuate errors 
in centering the charge vjithin the furance cavity and thus will 
prevent the establishment of axi-symmetric boundary conditions. 

Gradient Zone Length ; For any given temperature difference between 
the hot and cold zones, a decrease in the gradient zone length, X^, 
will increase the axial gradient, especially for large Bi (eq. 

[10]). h smaller Xg, hov/evor, will also produce larger radial gra- 
dients in the gradient zone,^"^^ and the precision with which the 
interface must be localized in order to satisfy radial gradient 
criteria will increase. 

Charge Diameter ; For a given system, any decrease in charge diameter 
will decrease Bi and thus the nondimensional axial gradient, but 
will increase the dimensional axial gradient (see Section 3.4). 

This effect is most pronounced for large Bi and when accompanied 
by a decrease in Xg. 

Furnace Temperature Difference Between Hot and Cold Zones ; Increas- 
ing will produce a proportional increase in the axial gradi- 

ent, eq. [11] . 

Crucible ; To prevent a severe reduction in the thermal con- 

ductivity of the crucible should be close to that of the charge. 

In systems with large conductivity differences between the liquid 
and solid charge, the conductivity of the crucible should be chosen 
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so that the nondimonsional axial gradient is maximized according 
to eg. [12], A thin crucible (5 close to unity) also serves to 
keep high. If vapor pressure considerations arc important, 

the resultant decrease in crucible strength can bo overcome by 
pressurizing the furnace system. An alternate approach to maxi- 
mizing consists of the use of coated metallic crucibles. 

Growth Rate ; If the generation of latent heat is significant (eg, 
[17]), a decrease in the growth rate (smaller P®sB|j) will serve to 
increase the axial gradient in the liquid at the interface. 

The present study of functional relationships concerning the 
axial tomporature distribution in Bridgman configurations is of 
obvious importance for gradient control and the related morpholog- 
ical stability of the growth interface. Moreover, axial gradient 
control is interrelated to the nucleation and segregation affecting 
radial thermal gradient control, dealt v;ith in Part XX of this work. 
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APPENDXX 


For an infinite charge length and crystal-melt interface posi 
tion within the gradient zone, the model of Section 2 yields the 
following expression for the axial gradient in the liquid at the 
interface ; 


m, 


de 

a? 
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+ POgRjj 
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= I’eclct number based on liquid properties 
PCg = Poolet number based on solid properties 
Bij^ « Biot number in hot zone (based on liquid properties) 
Bi^ = Biot number in cold zone (based on solid properties) 
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5 ■ dimonsionless interface position with respect to 

the center of the gradient zone, -y < < +y. 


The relation between interface temperature (0 ) and position 

m 

( 5 ^) is given by: 
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Note that an iterative solution is required if it is desired to 
determine given 0 . 

For a symmetric system, eg. [9] holds and m^ and m 2 can be 
approximated by eg. [20] . This leads to a simplified form of the 
axial gradient expression, which is used in Section 3. 



, qracUonts aU Uho crysUal-molt inUoirfaGo a 
lugh cqa, [10] and [24], dF/dc? and avorago 
, gradients (dO/dc)jj^, as determined from 
isional model of Pu and Wilcox. 
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FIGURE CAPTIONS 


Pig. 1 
Pig. 2 

t 

Pig. 3 


Pig. 4 


Pig. 5 


Pig. 6 


Pig. 7 


Pig . 8 


Pig. 9 


Fig, 10 


Modiiiod Bridgraan-Stockbarger configuration. 

Thermal model of furnace and charge with crucible. 

Effect Of thermal conductivity of a crucible on the 
effective Biot number. 

Effect of charge motion on the axial temperature dis- 
tribution in a charge. 

Effect of thermal coupling between furnace and charge 
on the axial temperature distribution within the charge. 


Effect of gradient zone length and Bi on the axial tem- 
perature gradient in the melt at the crystal-melt inter- 
face . 


Effect of unequal hot and cold zone Bi on the axial 
temperature distribution in the charge. 

Effee.t of charge diameter on the axial temperature 
gradient at the crysstal-roclt interface. (The axial 
coordinate is dimensional.) 

Effect of difference in thermal conductivity of crystal 
and molt on the axial temperature distribution in the 
charge. 


Effect of thermal conductivity ratio (R^^) and Bi on the 
axial tempai*aturo gradient in the melt at the crystal- 
melt interface. 
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Fig. 11 Effect of the generation of latent boat at the growth 
interface on the axial temperature distribution in the 
charge . 

I^ig, 12 Effect of charge position within the furnace on the 
axial tempeirature distribution. 
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